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Abstract: Radiation propagating over cosmo logical distances can probe 
light weakly interacting pseudoscalar (or scalar) particles. The existence of a 



spin-0 field changes the dynamical symmetries of electrodynamics. It predicts 
^ . spontaneous generation of polarization of electromagnetic waves due to mode 

mixing in the presence of background magnetic field. We illustrate this by cal- 
culations of propagation in a uniform medium, as well as in a slowly varying 
background medium, and finally with resonant mixing. Highly complicated 
correlations between different Stokes parameters are predicted depending on 
the parameter regimes. The polarization of propagating waves shows inter- 
esting and complex dependence on frequency, the distance of propagation, 
coupling constants, and parameters of the background medium such as the 
plasma density and the magnetic field strength. For the first time we study 
the resonant mixing of electromagnetic waves with the scalar field, which 
occurs when the background plasma frequency becomes equal to the mass of 
the scalar field at some point along the path. Dynamical effects are found 
to be considerably enhanced in this case. We also formulate the condition 
under which the adiabatic approximation can be used consistently, and find 
caveats about comparing different frequency regimes. 
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1 Introduction 



Light, weakly interacting spin-0 particles are predicted by many extensions 
of the Standard Model of Particle Physics. Such field might arise as pseudo- 
Goldstone bosons of some spontaneously broken chiral symmetry (PQ) [JJ, El 
[3] , in supergravity theories pfl [5] and the low energy string action P, [3, [8] 
as the Kalb-Ramond field. The mixing of light with pseudoscalar fields in 
propagation has a long history. Karl and Clark |9J looked for "cosmological 
birefringence" more than 20 years ago. Pseudoscalar-electromagnetic field 
mixing has long been explored by several authors [TQl [Til 131 HH EES EE] ■ 
Our focus is to explore use of polarization observables that can separate 
models. Our calculations find that polarization evolution can sometimes be 
highly complicated and not easily estimated by simple dimensional argu- 
ments. Placing limits based on rough arguments and dimensional analysis 
is far from reliable. Recently Csaki et al [T7] raised the possibility that su- 
pernova dimming might be caused by photon-pseudoscalar (7 — 0) mixing, 
re-igniting interest [18J in this problem. Whether or not supernovas dim, the 
overall physics of propagation turns out to be quite interesting, and the sig- 
nals of polarization particularly sensitive to careful treatment. The variety 
of effects is huge. As a result, one must question the validity of placing lim- 
its based simply on dimensional analysis or naive generalities. Instead, the 
logic has to be reversed to seek particular effects contradicting conventional 
physics wherever they occur, and then investigate the causes. 

There is another reason to consider polarization as a cosmological observ- 
able. Physical theories are not tested by fitting undetermined parameters: 
they are tested by testing their symmetries. The framework of general relativ- 
ity (GR) has certain symmetries by construction. Among these symmetries 
is duality, which predicts that the polarization of light in free space does not 
change. By the principle of Equivalence, GR then prescribes parallel trans- 
port of polarizations that are trivial in freely falling coordinates. It follows 
that observations inconsistent with duality might rule out perturbations of 
the metric as a false signal. Another symmetry of unitarity prevents the 
spontaneous appearance of polarization in an unpolarized beam of light. If 
such features contradicting GR are observed, then it rules out GR as an 
explanation of them, lifting a "degeneracy" of interpretation, and putting 
responsibility uniquely on something new. 

There is one spectacular signal. A spin-0 field generically induces po- 
larization even in a completely unpolarized beam. Indeed such a field can 
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spontaneously polarize the CMB itself. The effect is simple and quite distinc- 
tive. We are emboldened to suggest that the absence of polarization effects 
in propagation over cosmological distances, or in CMB derived quantities, 
would put stringent limits on the existence of a light spin-0 field. In this 
regard we have eagerly awaited the release of WMAP data on CMB polar- 
ization while preparing this article. We feel that the time for discovery with 
polarization observables may be ripe with this and other sources of data 
coming in the future. 



2 Symmetries of Light in the presence of spin- 
field 

The existence of a light scalar (or pseudoscalar) field changes the symmetries 
of the propagation of light. We consider a spin-0 field (ft coupled to the 
electromagnetic field strength F^ u by the action 



S = d A x\/—g 



-F^+g^e^F^F^ 
+g'^F^F flv + + - l -ml4> 2 + V{<t>) 



(1) 



We include a coupling to a current for completeness, while in practice this 
term will lead to modified propagation such as plasma frequencies, etc. For 
the purposes to be developed the potential V(4>) can be ignored as a small 
perturbation, and the metric g be replaced by a given background form. Note 
that we have included coupling of 4> to both e llua pF^ v F alS and F^F^, since 
both are allowed by the symmetries of the theory. If is a pseudoscalar, 
as we assume, then the coupling g'^ breaks parity symmetry. However this 
approximate feature of the low energy world is not a symmetry of Nature, 
being broken by the weak interactions physics. Even if the coupling g'^ is 
tuned to zero at some scale, a coupling will be induced by radiative correc- 
tions. In the rest of the paper, however, we ignore g'^ since the limit on this 
coupling is much more stringent compared to the limit on |19j . 

Free space electrodynamics has an interesting symmetry called duality. In 
the Hamiltonian density 7i in natural units, the electric and magnetic fields 
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E and B can be rotated: 

E' \ / cos# sin0 \ / £ 
5' y y -sin6> cos6> J \ B 

H = ^E 2 + ^B 2 (2) 

Although a symmetry of the Hamiltonian density, duality is not a symmetry 
of the Lagrangian density, which changes by a pure divergence. Noether's 
theorem represents this faithfully: the duality current is not conserved, 
but its divergence is an identity: 

K, = e^A"E al3 ; 
•9^ = e^ af3 F^F af3 . (3) 

As a consequence of duality symmetry, e^apF^F^ ~ E ■ B, which is not 
invariant, must be zero in a free-space plane wave. Otherwise symmetry 
would permit rotating it into a different magnitude, which is a contradiction. 
The duality current is deeply related to the helicity density: 

]-K = i-.A* k x E k , (4) 
z to 

where k is the wave number, u the angular frequency and Ax E is the canon- 
ical spin density^ Thus duality symmetry implies separate conservation of 
the spin and orbital angular momentum of light in a k eigenstate. Now the 
coupling of to e^apF^ 'F al3 breaks duality symmetry, allowing the helicity 
of photons to mix, and the plane of polarization of light to change in prop- 
agation - at least in the general case. Remarkably, the coupling of Eq. [1] 
vanishes for a single plane wave: an isolated photon is "safe" . However light 
cannot be protected under all possible conditions, and below we explore the 
effects of a background magnetic field. 

Let us contrast the 0-coupling with propagation in the usual cosmological 
model. In that event there is a coupling of light and gravity: the metric enters 
as an "optical medium" in the form of 



S = d^x\[^g 



--F^9 ua (x)9^(x)F aP 



(5) 



1 Strictly speaking, photon spin depends on the gauge choice. The interpretation comes 
from the gauge A = 0, which is nearly unique in having valid canonical spin and orbital 
angular momenta. 
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The energy of this system can be written 

H = ^(D-E + B-H); 

D k =g M {x)g» k {x)F, u] 
H k = \e ijk g^{x)g» k {x)F„ v . (6) 

This system has a new duality of rotating D and H into E and B. In a 
coordinate system where the metric is trivial, then D = E and H = B and 
the state of polarization is preserved. 

GR also maintains a symmetry under which the power transmitted scales 
with kinematic dependence on the scale factor, while the plane and magni- 
tudes of polarizations are preserved, in the sense of parallel transport along 
geodesies. The prediction of GR effects then consist of a coordinate trans- 
formation. In the isotropic, homogeneous metrics assumed in cosmology, the 
observable result is a red-shift. 

Free-space symmetries map into GR symmetries due to the principle of 
equivalence. Thus conclusions about them do not hinge on particular solu- 
tions. With GR serving as a "null" theory, the effects of 7 — <fi mixing can be 
separated and distinguished: simply observe the spontaneous polarization of 
light during propagation, for instance. 

2.1 Mixing Light with Pseudoscalars 

Astrophysics constraints limit the coupling to be very tiny, yet the enor- 
mous length scales of cosmology allows small cumulative effects to develop 
into large ones. Proceeding requires discussion of scales so that we can lin- 
earize Eq. CD to discuss mixing of modes in propagation. There are two 
important classes: background <ft, where the electromagnetic field is solved 
with a fixed pseudoscalar field, and background B, where a combination of 
the electromagnetic field and <fi are solved for a given magnetic field. Here 
we set the background (f> field to zero and put our effort on the second case. 

We will use symbol B for the magnetic background, and B for the field 
in propagation: the total magnetic field is B + B. Although galactic and 
intergalactic B ~ fiG fields are small, their effects are likely to be large 
compared to the effects of <fi. A quick calculation shows that E/c and B 
of a typically weak but observable signal, such as light from a QSO or the 
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CMB, are also small compared to pG. We seek to linearize consistently with 
E/c « B, and B « B. There is a well-defined 3-state propagation problem 
of mixing 2 light polarizations with coupling to <fi. We first obtain the non- 
covariant form of Maxwell equations, as follows, and with no approximations: 

V-E = g^(P-(B + B)+p- (7) 
Vxi + ^ = o; (8) 

VxB-^ = g^ExV^-iB + B^+j; (9) 

V-(B + B) = 0. (10) 

Here £>« + Bi = \tijkF^ k and E\ = F° l are the usual electric and magnetic 
fields. Besides this we have the equation for the pseudoscalar field 



<9 2 



- V 2 <P + mi<P = -g^E-(B + B) (11) 



dt 2 

Gauge invariance is explicit and one can check current conservation directly, 

The dynamics of this coupled 3-field system are complicated and been visited 
by various approximations in the literature. 

We assume B solves the zeroeth order Maxwell equations with no <fi back- 
ground. Our linearized equations for E/c « B, B « B are 

V-E = g^-B + p; (12) 
OB 

VxE + — = 0; (13) 

V-B = 0. (15) 

These equations remain exactly consistent with current conservation, V • j + 
dp/dt = 0. 

Proceed to get a wave equation for E by taking the curl of Faraday's Law, 
and substituting into Eqs. IS Replacing B ~ B gives 

- V 2 ^ + — = g+g^Z - ^V(V0 ■ B). (16) 
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— » — » — * — * 

In this equation the longitudinal part of E mixes with V0 • B. That is, E is 
not perfectly transverse in general when coupled to 0. 

Continuing along these lines towards a general solution greatly increases 
the complexity of the equations. We will be content here to show illustrative 

— * — » — » 

calculations restricted to the limit V0-£> = 0. In this case E is transverse and 
the system readily collapses to 2-state mixing. We have extensive calculations 
for the more general cases, but we feel this paper would not be improved 
by introducing long and complex calculations. Our detailed calculations 
show that the numerical importance of Bl ^ effects for the parameter 
regions considered in this paper are negligible. The only exception to this 
may possibly be in the case of resonance where the correction terms may be 

— * — * 

lumped together with other parameters in any event. The virtue of V0-£> = 0, 
even if it is not the most general case in practice, is that one can get a feel 
for the rich interplay of several dimensionful scales including g^, B, m?, the 
plasma frequency u> p , and the propagation distance z. The amazing range of 
phenomena one should seek observationally is well illustrated. 

Requirement 0: Existing Physics Known physics of electromagnetic 
propagation in matter must be taken into account. It is well known that the 
plane of linearly polarized light rotates in a propagation through magnetized 
plasma: the Faraday effect. Our symmetry arguments are consistent, because 
duality is broken by the coupling to electric (as opposed to magnetic) charges 
and currents. Fortunately Faraday rotation is quite frequency dependent, 
and most important for low frequencies below the GHz regime. 

We also take into account the practical "photon mass" in the form of 
the plasma frequency. It is not an effect that can be ignored, and we do 
incorporate it in the momentum-space propagation equations below. 

Requirement 1: Phase Accuracy Observables depend on the wave 
number differences Ak accumulating phase AO by propagating over large 
distances Az. An accurate approximation needs the absolute phase error 
(symbol 5) to be small compared to n, namely 5(Ak±Az) « 1. This 
requirement becomes more and more demanding as Az — > oo. Everyone 
recognizes this requirement, listed here for completeness. In the absence of 
degeneracies, the wave numbers of propagating waves can often be obtained 
consistently in powers of g^, provided the corrections are controlled: see 
below. 
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Requirement 2: Gauge Invariance It is very important for all symme- 
tries to be respected, and in particular, any violation of gauge invariance is 
unacceptable. Fortunately our equations set up in a gauge invariant manner 
coincide with the standard literature in the limit of B ■ B = 0. 

Requirement 3: Respect for Limit Interchange Shortly we will solve 
the propagation and present a common series expansion in g| to simplify 
the formulas. We will find that the series contains factors of the frequency 
uj. The limit of g^ fixed and small, and taking u — > oo, does not commute 
with uj fixed and g^ — > 0. Since the frequency of a eV photon in units of the 
cosmological length scale is huge, while the coupling constants contemplated 
for is tiny, one needs to examine series expansions carefully to be sure they 
apply to the problem being solved. 

3 Pseudoscalar-Photon Mixing: Uniform Back- 
ground 

Choose the coordinate system with the z-axis along the direction of prop- 
agation of the wave and the x-axis parallel to the transverse component 

— * 

of the background magnetic field Br- Seek solutions with harmonic time 
dependence e~ lwt . Denote the component mixing by subscript || and its per- 

— * — * 

pendicular complement by subscript _L. Define symbol A = E/to, which 
simplifies equations much like the vector potential, except there is no "i" 
and A is gauge invariant. 

The perpendicular component Aj_ does not mix with <fi. The wave equa- 
tion for the mixing of A\\ and (f> can be written as, 

(» 2 + ®{t))- M {t))=° (17) 

where the "mass matrix" or "mixing matrix" is 




and Br is the magnitude of the vector Br- At this point we took into account 
the plasma frequency u^, which is a (non-local) gauge invariant mass for E. 
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Transform to a new basis 



<K*) J V 4>{z) 



where is the orthogonal matrix 

O = 



cos 9 — sin 9 
sin 9 cos 9 

which diagonalizes the mixing matrix in Eq. [TTJ Diagonalization gives 

2g< t> uB r 



tan 29 - 



with mass eigenvalues 



We list the leading order expansion in g^, by which 

2 _ 2 g Br u 

2 2 

M- = 2 2 • 

We noted earlier that examination of scales is needed to justify this 
Consistency requires 

gBru <<L 



Continuing, the equations are written in the diagonal basis as 

(uj 2 + 9 Z 2 )Z|| - = o 

(U 2 + d 2 z )(j) - /i^ = 

which can be easily solved to give 

A\\(z) = Z||(0)e 4{w+AA)2 
$( z ) = 0(O)e^ +A ^ )z 



and 

,2 _ ,,2 ~ - 



w + A A = \/uj 2 - ^ « uj - ^± 



uj + A^=^uj 2 - f i 2 _^uj~^- (28) 

We will always work in the limit uj 2 » justifying retention of just the 
leading terms. The phase difference A^ — A^ in this limit is found to be 



A * " Aa * 2u~j V (tt * " m ' )2 + ^Brujf (29) 
The perpendicular component A± does not mix with and is given by 

A±{z) = A ± {0)e ikoz (30) 
where ko = sju) 2 — uj 2 ~ uj — uj 2 /2uj. 

3.0.1 Typical Scales and Units 

We use parameters typical of the Virgo supercluster for illustration. Gener- 
ally the plasma frequency is given by 

ujI = = - fi e - 3.7 x 10~ 14 eV) . 31 

p m e 10- 6 cm- 3 v ' v ' 

For the supercluster, the plasma density n e ~ 10~ 6 cm" 3 and the magnetic 
field is of the order of 1 B ~ 0.1 /xG coherent over a distance scale of order 
10 Mpc [20J. As mentioned earlier, use of typical B values in Br formulas is 
done to take advantage of the transparent simplicity of that limit. 
In the limit of « uj p , there is a useful length scale I defined by 

1 = ( 2^ « in6 ptt °- 04 M P C ( 32 ) 

(uj 2 — mi) 10 b GHz 

where u — j-. We refer to this as the oscillation length in analogy with a 
similar variable in neutrino physics [2T] . The value 10 6 GHz ~ 4 eV is a 
handy upper order of magnitude for optical frequencies. 

A typical upper limit [22], [23] EU [25] on the coupling parameter is 
6 x 10~ n GeV" 1 . In Fig. [T] below, we show results for g$ ~ 10~ 12 GeV -1 , 
more than one order of magnitude smaller than the limit. For a magnetic 
field of 0.1 /iG we find that the product g^B can be expressed as g^B = 0.215 
Mpc" 1 . For convenience the scales are summarized in different units in Table 

m 
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quantity 


typical value 


alternate units 


B 


0.1 ii G 


4.78 x 10 49 Mpc~ 2 


9d> 


10- 11 GeV- 1 


6.4 x 10" 50 Mpc 


u; p 


3.7 x l(T 23 GeV*/ in e e 3 

V 10 °cm 3 




OJ 


10~ 5 - 1 eV 


1.6 xlO 24 - 1.6 x 10 29 Mpc" 1 



Table 1: Typical values of dimensionful scales in different units. If not other- 
wise specified we use ft, = c = 1. The value of listed is far below published 
limits of < 6 x 10~ n GeV" 1 [22]. 

3.0.2 Small Mixing Limit: a Simple Example 

For an example we examine the most innocuous case of mixing angle 9 — > 0. 
Specifically, this is the limit 

m 

T «1. 



\m 2 , — ojI 



By Eq. [251 the limit of small mixing is just the same limit in which the naive 
Taylor expansion of Eq. [24] in small ~ applies. In this limit Eu, E± 
and are the approximate propagation eigenstates. One might think there 
are no observable effects, and in particular, there would be no substantial 
"dimming" of intensity. 

However there is an important relative phase shift in propagation: 

E||(^) = E||(0)e^ +A ^; 
E±(z) = E ± (0)e ikoz . 

The physically observable density matrix p is given by 

p ~ V < E ± E* > < E ± E* ± > ) ' (33) 

where < > denotes the statistical averages occurring in propagation. Under 
coherent conditions this result predicts a cumulative rotation of the plane of 
a linear polarization due to the off-diagonal term: 

E\\El{z) = £||£l(0)e i( ^ /2 ^+ /2w)2 . (34) 

So long as the mixing is small, we may insert the g 2 expansions of Eq. [241 and 
predict the angle of rotation increases linearly with frequency. Yet there is 
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always a limit in which the frequency is large enough to cause strong mixing. 
We explore this for m| « u^. 

Let the magnitude of the phase angle between E\\ and E± be denoted \- 
From Eq. [34] we have 



Using g^Br ~ Mpc -1 as a typical value, the phenomenon of a rotating po- 
larization might be readily observed to be linear in u at radio frequencies for 
decades above the GHz regime. 

Approaching the optical region, both the formula for small mixing and 
the expansion of Eq. [24] break down, leading to very interesting possibilities. 

3.0.3 Comments on Propagation 

It is interesting to contrast the results above with propagation in a disper- 
sionless and non-dissipative medium, such as "free-space" with gravitational 
fields. It is obvious that the intensity of a wave (Stokes I) is preserved up 
to the kinematic redshift of propagation 0. It is less obvious that the degree 
of polarization (Stokes Q/T) cannot be changed by any 2-state purely elec- 
tromagnetic propagation in a dissipationless medium. The origin of this is 
"unitarity" of propagation which can be written as a unitary evolution oper- 
ator, just as in quantum mechanics [26]. As a consequence of 2-state electro- 
magnetic unitarity, a particular circular polarization (say) can be converted 
into a linear one, but a linear polarization cannot be made to come from an 
unpolarized ensemble. In the weak mixing region, an unpolarized density 
matrix (a multiple of the identity) evolves with no change whatsoever. This 
means that basing observations or parameter limits on unpolarized quantities 
can easily miss effects that polarized observables would readily detect. 

Spontaneous polarization is more than rare, and in pure electrodynamics 
any form of polarization is usually associated with a corresponding extinc- 
tion. An ordinary polarizer plate is typical, with the degree of polarization 
scaling directly with the degree of extinction. The same goes for common 
sources of polarization such as Compton scattering invoked in astrophysics. 

2 Reflections, namely the generation of backwards moving waves, are always possible in 
any varying medium. We work in the limit that they are negligible. 



X 




X < 



zg^Br- 



(35) 
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Momentarily we will discuss stronger mixing cases where mixing of light with 
pseudoscalars can lead to spontaneous polarization. 



+ 



sin 29 cos 29 — sin 29 cos 29 coslzfA^ — 



3.1 General Density Matrix 

We can now determine the general density (coherency) matrix elements after 
propagation through distance z in terms of the matrix elements at the source. 
Inasmuch as we can interchange A with E/ou, it is convenient to report density 
matrices of A. Then: 

< A\(z)A\\{z) > = \ < Af(0)A||(0) > [l + cos 2 2fl + sin 2 2flcosKA^- A A 

+ - ^ (0 ^ (0) > [ S in 2 29 - sin 2 29 cos[z(A^ - A A )}} 

^ <A|(O)0(O)> 
I 2 
— i sin 20 sin [^(A^ — A^)] +c.c.| 

We gave the most general expression above keeping all the correlators at 
z = 0. One might think it reasonable to assume that < A*^ (0)0(0) > is zero. 
However in considering propagation through an intergalactic medium one 
expects a large number of magnetic domains uncorrelated with one another. 
The general expression is needed to describe sequentially the propagation 
through different domains. Given the randomness of such processes, the 
generic situation is one where all the correlators on the right hand side will 
be nonzero. 

We also find, 

< A\{z)A ± {z) > = (cos 2 9e iFz + sin 2 9e iGz ) < A|(0)A ± (0) > 

+ cos0sin0(e iF2 - e lGz ) < 0*(O)A ± (O) > (37) 

and < A* ± (z)A ± (z) >=< A* ± (0)A ± (0) >. The phase factors F and G are 
given by 

F = u 



G = uj 







A- 










A- 


V- 

uo 2 


1 , 

m — (u_ — to 
2u 
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Here again it is reasonable to expand F and G and keep only the leading 
order terms in u 2 /uj 2 and fi±/u 2 . Using these correlation functions we can 
compute reduced Stokes parameters I,Q,U,V (or So, Si, S2, S3) 

I = < Af(»A||0) > + < A* ± (z)A ± (z) > 

Q = < A\(z)A\\(z)> - < Al(z)A L (z)> 

U = < A\(z)A L (z) > + < A* ± (z)A^z) > 

V = i{- < A\{z)A±(z) > + < A* ± (z)A { \(z) > (38) 

at any position z. The standard Stokes parameters of the same name are 
simply u 2 times the above. We will often remove the scale of intensity, 
normalizing I = 2, as in an unpolarized wave |^4||| 2 = |^4_i_| 2 = 1- 

3.1.1 Special Cases 1: Spontaneous Appearance of Polarization 

Suppose the initial beam is unpolarized and the initial correlator < At (0)0(0) > 
is also zero. We scale out normalizations so that 

< Af(0)yl||(0) >=< A*j_(0)A±(0) >= 1 

and assume that all other correlators vanish at z — 0. The expressions 
simplify to: 

< A\{z)A\\{z) > = i (l + cos 2 2# + cos[,2(A (/ , - A A )]sm 2 29) + ^ 

x < 0*(O)0(O) > (1 - cos[2(A - A A )\) sin 2 29 (39) 

< Al(z)Aj_(z) > = (40) 

< A* ± (z)A ± (z) > = 1 (41) 

This wave spontaneously acquires a linear polarization oriented along B 
during propagation. No circular polarization is developed. We plot the 
degree of polarization in Fig. [T]for initial conditions < 0* (0)0(0) >= 0. The 
degree of polarization p here is given by 

JQI \ < Al(z)A\\(z) > - < A* ± (z)A ± {z) > \ 
P I < Al(z)A\\{z) > + < A* ± (z)A ± (z) > ' ( j 

given Stokes U = V = 0. The degree of polarization accumulates to a sizeable 
magnitude and could produce observable consequences even for exceedingly 
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Figure 1: The degree of polarization p as a function of the length parameter 
/ (Eq. 1321) for uniform background propagation for two different value of the 
coupling parameter g$. Values of are given in units of GeV -1 ; B = 0.1 \i 
G, n e = 10 -6 cm -3 . The length scale I is in units of Mpc (z//10 6 GHz). The 
propagation distance is taken to be 10 Mpc. 
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small couplings. In Fig. [T]the degree of polarization is shown for two different 
values of the coupling parameter = 1 x 10~ 12 , 2 x 10~ 12 GeV -1 . The 
background medium parameters are chosen to be B = 0.1 fiG and n e = 10 -6 
cm -3 . The propagation distance is taken to be 10 Mpc. This is also the 
distance over which the background magnetic field is assumed to be coherent. 

The experimental signature of this mixing for a uniform background 
would require fitting the observed degree of polarization to a number of 
parameters. As shown by the previous examples, the extrapolation between 
radio frequency and optical frequency observations needs care and attention 
to series expansions. Nevertheless the spontaneous polarization due to mix- 
ing is clearly distinguishable from extinction, which is expected to show a 
simple monotonic increase in the degree of polarization with frequency. 

3.1.2 Special Cases 2: Initial Polarization 

Consider now the case where the incident beam is partially polarized. A 
sample of results for randomly chosen intial conditions is shown in Fig. [2j 
The Stokes Q/I, U/I and V/I parameters are presented as a function of the 
length parameter I, equivalent to the frequency of the incident wave. Initial 
conditions are < 0*(O)0(O) >= 0, < Aj(O)0(O) >= and < A* ± (0)<p(0) >= 0. 
A more general result is discussed later. One interesting pattern that emerges 
from this figure is that U and V show a much larger variation with frequency 
compared to to Q. This is related to observations of Ref. [27] that at low 
frequencies the dominant physical effect arises due to the phase in A\\, an evo- 
lution of the polarization direction of the wave, which is permitted by 2-state 
unitarity. The figure supports near conservation of \JU 2 + V 2 which is rela- 
tively independent of frequency, both in the low and high frequency regime. 
Indeed this variable shows spectral dependence only in a relatively narrow 
region. In the high frequency regime, of course, all the Stokes parameters 
show dependence on frequency. 

A possible signal of the pseudoscalar mixing is the correlation between 
the different Stokes parameters, in particular between U and V. In Fig. [3] we 
show a sample of results to illustrate the correlation between the Stokes pa- 
rameters U and V for some randomly chosen parameters and initial state of 
polarization. We find these Stokes parameters do show an observable depen- 
dence on each other. The dependence of U onV is found to approximately 
follow an ellipse for large frequencies. If the frequency is small such that 
the oscillation length I << L, where L is the propagation distance, then 
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Figure 2: Reduced Stokes parameter ratios Q/I, U/I and V/I as a function 
of the length parameter I in units of Mpc, (Eq. [521 . for uniform back- 
ground propagation with = 2 x 10~ 12 GeV -1 . The degree of polarization 
p = \/Q 2 + U 2 + V 2 / 1 and the variable y/U 2 ~+~V 2 /I are also shown. Other 
parameters are B — 0.1 p, G, n e = 10 _6 cm -3 , Propagation distance =10 Mpc, 
m^/ujp = 0.1. The initial state of the polarization has been chosen randomly 
such that Q/I = 0, C//J = 0.4 and V/I = 0.1. 
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we find that the dependence of U on V is considerably more complicated. 
Remarkably, the Stokes parameter Q (degree of polazarization) also shows a 
simple dependence on U and V (polarization direction controllers) at large 
frequencies. This is shown for some randomly chosen parameters in Fig. |H 
These illustrative calculations show that a fascinating range of physical 
phenomena are generated by the mixing of light and pseudoscalars. 

4 Pseudoscalar-Photon Mixing: Slowly Vary- 
ing Background 

For the rest of the paper we will consider cases in where the background is 
slowly varying. There exists a large literature regarding the integration of 
differential equations with slowly varying parameters. The reader interested 
in the mathematical problem can consult the classic reference of Olver [28]. It 
is an understatement to say that the variety of situations is huge. Rather than 
make sweeping statements, we will use a workmanlike series of definitions and 
approximations that we check are suited to the limits we study. 

First restrict attention to the simple case where the direction of B is uni- 
form. Assume that the background varies sufficiently slowly to justify neglect 
of terms involving the derivatives of the background plasma density (or the 
magnitude of the magnetic field). We thus enter into the adiabatic approxi- 
mation, a general branch of mechanics, which holds that the amplitude in an 
eigenstate is relatively constant when the parameters defining the eigenstate 
change sufficiently slowly. We again obtain two decoupled equations Eq. [26] 
which can be solved by using the ansatz 

A\\{z) = ~A ll (0)e iujz+l ^ AAdz ' 
$(z) = 0(O)e 4W2+ ^o A ^ 2 ' (43) 

The phase factors are given in Eq. [281 

In the case there is a conservation law - the conservation of transmitted 
energy flux, similar to the conservation of probability in quantum mechanics - 
the WKB approximation would indicate certain pre-factors of order l/y/k(z), 
where k(z) are the exponential integrands. We review this with an iterative 
WKB procedure developed in Section 5 which controls coefficients in a certain 
expansion. However extensive calculations support wide applicability of the 
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Figure 3: A sample of results showing the correlation between the Stokes 
parameters U and V for randomly chosen parameters and initial state of 
polarization under propagation in a uniform background. Parameters (in 
arbitrary units) are (a) g^B = 2, L = 10, 0.04 < I < 20, (b) g^B = 
10, L = 10, 0.4 < I < 300, (c) g^B = 1, L = 50, 0.2 <l < 800 and (d) 
g^B =10, L — 10, 0.04 < / < 100; mj/wj = 0.1 for all the plots. A simple 
correlation is seen only for frequencies larger than a minimum frequency. At 
low frequencies the relationship becomes very complicated. 
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Figure 4: A sample of results showing the correlation between the Stokes 
parameters Q and U. Parameters and initial states of polarization are taken 
to be same as used in Fig. [31 The oscillation length I is (a) 1.1 < I < 20, (b) 
0.4 </ < 300, (c) 1.7 < I < 800 and (d) 0.27 < I < 100. The lower limits 
on / are taken slightly larger than the values used in Fig. [3] in order that a 
simple relationship emerges between Q and U. 
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simpler formulas, inasmuch as k(z) tends to cancel out of polarization ratios 
and not to develop cumulative run-outs in many circumstances. 

Using these results we can now evaluate the electromagnetic field at any 
position. For this purpose we express the fields A(0) and 0(0) in terms of 
A||(0) and 0(0) and the mixing angle 0(0). The final expressions for the fields 
A\\(z) and <f>(z) are given by 



A\\(z) = e l "M||(0) 

+ e^ 2 0(O) 

(f)(z) = e^ 2 A||(0) 

+ e^ 2 0(O) 



cos 9 cos 9 exp yi J A^dz' J + sin 9 sin 9 exp yi J A^dz 
cos 9 sin 9 exp J AAdz'^j — sin 9 cos 9 ex P J A^dz' 
sin 9 cos 9 exp J A^dz^j — cos 9 sin 9 exp J A^dz' 



sin 9 sin 9 exp ( i \ A^dz' j + cos 9 cos 9 exp J A^dz' 



(44) 



where on the right hand side 9 = 9(z) and 9 = 9(0). The resulting expres- 
sions for the coherency matrix elements is given by 

< A|(0)Aii(0) > 

< A\(z)A\\(z) > = " [1 + cos 2g cos 2g + sin 2fl sin 2fl cos $] 

+ - - > [1 - cos 2# cos 2^o - sin 29 sin 2# cos $] 



+ 



4(0)0(0) > 



+ i sin 26> sin $ 



+ c 



cos 2$ sin 29 — sin 26 1 cos 29 cos $ 
.c.) (45) 



where 



$ = / (A^ — A^)dz' 
Jo 



(46) 



< A^A^) > = < A*(0)A ± (0) > cos cos 6> exp ( — / dz'\p+-u). 



+ sin sin exp ^ dz' [/ii - u;^ j 
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+ < <p*(Q)A ± (0) > cos0sin0 o exp Qjjf dd\pL\-u§\ 

— sin cos $0 ex P {^2~ J dz'\p?_ — • (47) 

The above result is valid as long as the medium changes slowly so that the 
term 

*>*(«.<»«. = "o 1 )^) («) 

is negligible compared to mass term 

>i \ /A* 



/i 2 _ 

We are justified in ignoring the derivative of the mixing angle 9 if 

|*1 « § (49) 

The condition for adiabaticity will be formulated in greater detail in section 
5. There we estimate the transition probability from one local eigenstate 
to another. As we shall see, the condition given above, Eq. H9j is sufficient 
provided the difference between /i + and /i_ is of the order of these eigenvalues. 
If there is a delicate cancellation between these two eigenvalues then on the 
right hand side in Eq. H9] more care and specialized techniques may be 
needed. 



4.1 Resonance 

We next discuss the interesting case where the plasma frequency uj p becomes 
equal to the pseudoscalar mass somewhere along the path. Given the wide 
range of plasma frequencies observed in the astrophysics there may be many 
situations where such resonance might occur. Examples include propagation 
of light or pseudoscalars from a dense medium, such as an AGN, GRB or 
cluster of galaxies into the intergalactic medium. Another example includes 
propagation of light from pulsar magnetosphere into the interstellar medium. 
Let the radiation initially propagate from the region z = containing high 
plasma density u p > towards increasing z where a resonance occurs. 
Assume that initially the mixing angle 9 « 1, and recall 

tan20 = %^«l. (50) 

rr, 
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Let the plasma density decrease slowly along the path such that at the ob- 
servation point > uj p . Somewhere along the path at z = z* the right 
hand side of Eq. [50] becomes infinite, which we will implement with angle 
20(z*) -> tt/2 . To finally fix the initial conditions, let < Aj(0)A||(0) >^ 0, 
< A|(0)v4_l(0) >^ and < A* ± (0)A ± (0) >^ , that is, a generically mixed 
wave. The remaining correlators, which involve 0, are taken to be zero at 
z — 0. We compare the conditions after traveling a distance Az = L to the 
observation point. From Eq. [45],H7]we find: 

< A\{L)A\\{L) > « 0, 

< AJ(L)A ± (L) > ^ 0, 

< A* ± (L)A ± (L) > = < A* ± (0)A ± (0) > (51) 

if | tan 28\ << 1 at z = L. These are predicted independent of frequency and 
the distance travelled, provided the wave crosses the region where = u p . 
Almost independent of the state of polarization at the origin, the wave be- 
comes completely linearly polarized at the observation point. The orientation 
of the electric field vector is perpendicular to the background magnetic field. 



Studying resonance using a slowly varying approximation requires some 
care. The consistent requirement is that the assumption of slowly varying 
background be satisfied. We now make this assumption more precise in 
order to obtain a constraint on the derivative of the background plasma 
frequency. For this purpose we assume that the background plasma density 
varies linearly with position such that 

<-"l = C-az (52) 

where C and a are constants, z is the distance of propagation and u is some 
chosen value of the frequency of the wave. It is convenient to work with the 
rescaled variables, 

~2 U l ~2 m l ~ & 

u p = — ,m = — , u = — , (53) 

O>o O>o Co>o 



3 Generally all resonant systems have a special phase at resonance. We acknowledge R. 
Buniy for discussions and work long ago on similar resonant phase evolution in classical 
mechanics. 
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since for cosmological applications tu 2 has a value of order unity in units of 



Mpc -1 . The origin of the wave corresponds to z = and the observation 
point z = L. 

If at z = the plasma frequency is larger than pseudoscalar mass then 
the constants C and a are taken to be negative. We assume that the wave 
crosses the resonance region where — u p . As discussed earlier, in this 
case our approximation is valid if 

\ff\ « ■ (54) 

This is discussed in greater detail in the next section. This constraint reduces 
to 

\a\ « A{ g<i> Bfu (55) 

at the resonance point = u p (z). Now we come to the point of these num- 
bers: for the coupling — 6 X 1CT 11 GeV -1 the constraint is satisfied for 
the magnetic field, plasma density and the distance scale corresponding to 
the Virgo supercluster, provided the frequency of the wave to » 1.7 x 10~ 4 
eV. Since this is well below optical frequencies and well above the regime of 
GHz radio frequency studies, one should be prepared to see different phe- 
nomena in the optical and radio regimes. One can both use a slowly varying 
approximation, and observe resonant mixing consistently. 

Let us turn to the phases that appear in Eq. H7J for the case of linearly 
varying background. We find 

r dz' ^ - Up = \[-az 2 + 2zC) ± ^JAgmW + C 2 
Jo 4 4a 

± aZ ~ a ° ^jAg 2 B 2 Cu 2 + (az - C) 2 



( g<t) Buj 



~,\2 



a 



log ( az - C + \/(2g^BCu) 2 + (az - C) 2 

(56) 



± 



- log! -C+ ^(2g^Bu) 2 + C 

Since these phases are equal to f Q dz'[^\ — u 2 ]/2u we find that for a wide 
range of parameter space these phases vary roughly as 1/uj. The correspond- 
ing polarization observables change roughly inversely proportional to u. 

We can understand this behaviour of the phases appearing in the corre- 
lator < A*,(z)A±(z) > as follows. Consider the phase f Q dz\ji 2 + — uj 2 )/2uj, 
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which is basically the integral of the difference of the squared masses of one 
of the eigenmodes and the perpendicular component divided by 2uj, as ex- 
pected for relativistic particles. In the limit of small mixing the difference of 
the squared masses is proportional to the mixing angle squared and hence is 
proportional to uj 2 . In this limit the phase increases linearly with uj. How- 
ever if uj 2 — mi — > along the path there is a limit interchange. We need a 
different series expansion, for example 

A ^^ + l H+<) + ^^. (57) 

This contradicts perturbation theory in g^, because resonance is always a 
non-perturbative phenomenon, as signaled by the appearance of l/g^ in the 
expansion. In any event, given that occurs, the linear increase of phase 

flattens out near uj 2 — m| ~ 0, and can even decrease for large uj. 

So far we made the resonance occur. In the limit of small frequencies 
the phase factors become very large. For example in Eq. [56] the first factor 
contributes the term miz /2uj to the phase. For large propagation distance 
z and small uj this factor is clearly much greater than unity. Hence one 
expects that in general the Stokes parameters will show rapid fluctuations 
as a function of uj in this limit. However since g^B occur in a product, the 
small limit can again be reversed if E becomes very large! 

Continuing, we evaluate the phase $ appearing in Eq. H5] for the case 
of a linearly varying background plasma frequency, Eq. [52] Results can be 
extracted from the integrals given in Eq. [56] by using 

* = f\^-^W = [ g-gK (58) 

In fig. [5] and [6] we show two examples of the result expected in the case 
of resonant mixing. Here the initial state of polarization has been chosen 
randomly. The parameters values (in arbitrary units) given in the figure 
may be translated to those relevant for cosmological propagation by taking 
the units in the appropriate powers of Mpc. The pseudoscalar mass is chosen 
such that m 2 ^ = m^/ujo = 1, where ujq is defined in Eq. [52] As expected, 
observable effects are large at low frequencies compared to the corresponding 
results of a uniform background. We also find rapid fluctuations in this 
region. The analytic result, obtained in the adiabatic limit, is found to be 
in good agreement with the numerical result as long as Eq. is satisfied. 



25 




Figure 5: The Stokes parameters Q/I, U/I and V/I as a function of the 
length parameter 2u/(Up(z) — mV) (in arbitrary units) for the case of reso- 
nant mixing. Solid and dashed lines represent the results obtained by direct 
numerical integration and by using the analytic expressions Eqs. I45landl471in 
the adiabatic approximation. The analytic result is in very good agreement 
with numerical result as long as Eq. satisfied. Here &p(z) = m^—C+az; 
g^B = 0.1; mi = 1. Initial conditions at z = are p = 1.0, Q = 0.0, U = 1.0 
and V = 0.0. 
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Figure 6: Another example of Stokes parameters for resonant mixing. Pa- 
rameters g^B and L are approximately the same as Fig. [2] and = 1. 
Initial conditions are p = 1.0, Q = 0.0, C/ = 0.894 and V = -0.447. Solid 
and dashed lines represent the results obtained by direct numerical integra- 
tion and by using the analytic expressions Eqs. 05] and H7] in the adiabatic 
approximation. The analytic result is in very good agreement with numerical 
result as long as Eq. [55] is satisfied. 
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At smaller frequencies the two results start to disagree. As expected, we 
numerically find that the Stokes parameters approach their initial values in 
the limit u — *■ 0. 

In fig. [7] we show the dependence of the degree of polarization on the 
distance of propagation for the case of resonant mixing. Here the wave is 
taken to be unpolarized at source. The remaining parameters are taken to 
be g^B = 1.0, — 1, C — —0.5, a = —0.1 and the oscillation length at the 
observation point / = 2u/(Up(L) — m|) = 0.2. The figure clearly shows the 
large change in the degree of polarization as the wave crosses the resonant 
point. For comparison we also show the case of uniform background with 
essentially the same parameters as for the resonant case. The only difference 
is that the (uniform) oscillation length parameter I = 2uj/{uj'i — m|) = 0.2. In 
this case the wave acquires a very small degree of polarization in comparison 
to the resonant case. 

5 Adiabatic Approximation and Corrections 

In this Section we analyze the pseudoscalar-photon mixing for slowly varying 
background more systematically. This will give us an estimate of corrections 
to our results above. We assume the adiabatic limit applies: for any changes 
of parameters, we take derivatives such as dB/dz — > 0, with the propagation 
distance L —>■ oo, and the cumulative change AB fixed. 

Borrowing standard notation from quantum mechanics, we rewrite the 
basic equation, Eq. [T7J as follows 



We denote the eigenvectors of the mixing matrix as |n >, where in the case 
of resonant mixing of An with <f>, n = 1, 2. The eigenvalue equation can be 
written as 



Here the eigenvectors also vary slowly with z. The general solution to Eq. 
[591 can be written as 




(59) 



M\n >= fi 2 n \n > 



(60) 




(61) 



n 
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Figure 7: Degree of polarization as a function of distance for the case of reso- 
nant mixing. The wave is taken to be unpolarized at source. The parameters: 
g<t,B = 1.0, ml = 1, C = -0.5, a = -0.1, I = 2u/{u 2 p {L) - mj) = 0.2. The 
curve for uniform background (nonresonant mixing) is shown for comparison. 
All the parameters in this case are same as above, with the (uniform) mixing 
length I = 2uj/(uj 2 p - mj) = 0.2. 
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2 

where u n ~ u — ^ Substituting this into Eq. [59] an d taking the overlap of 
the resulting equation with < m\ we find 

d z a m = -l—a m - -iVW^^"^ < m|d> > (62) 

n 

We have dropped all terms involving two derivatives of the slowly varying 
quantities such as a n (z), \n> and the exponent. Using the eigenvalue equa- 
tion, Eq. [60] we find 

< m\(d z M)\n > 
<m\d z \n>= lK 2 J 2 l m^n (63) 

The resulting solution for be written as 

a m (z) = e-^^ dz ' {d ^ m)/uJm b m {z) (64) 

where b m (z) satisfies 



Urn ^ Hi,- [1% V u n{Z)U m (0) 

n,n^m V 

We pause to assess the result if we simply assumed b m (z) were constant. 
The pre-factor in Eq. [64] is integrated via 



exp (-1 J Q dz\d zr u m )/u m ) = J (66) 



which is the well-known WKB 1 / \Jk{z) factor cited earlier. In effect we have 
derived this factor and pushed the remaining derivatives in the equation for 
b m (z), Eq. [65] It is clear that this exponential term is approximately equal 
to one. 

We note that the equation for b m is very similar for the equation [62] for 
a n . This sets up an iterative scheme. Since slowly varying oj n ~ u, we can 
approximate Eq. [651 as 

d , bm „ yj b <™\fM)\-> e ,U^,,-«.,,, , (67) 
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plus small corrections in a series expansion of \uj — uo n \/uo « 1. For the 
cosmological applications, we consider here, the parameters ranges are such 
that the correction terms are extremely small. 

The value of b m (L) — b m (0) gives an estimate of the transition amplitude 
between the two local eigenstates after the wave has propagated a distance L. 
Let us now estimate this with dimensional arguments. In the adiabatic limit, 
b n , as well as the term < m\(d z M)\n > /(/^ — jtt£), varies slowly compared 
to the exponential. Hence we may take these to be approximately constant 
along the path. We can now integrate Eq. [67] to obtain the transition ampli- 
tude. We find that b m (L) — 6 m (0) is of order u/(IA/i 2 ) as long as we ignore 
the case of resonant mixing, where the two eigenmodes come very close to 
one another at some point along the path. Hence we find that in the limit 
L » u/Afi 2 the transition probability is negligible: this is the usual jus- 
tification for the adiabatic approximation when the internal dynamical time 
scale is very short compared to the time scale for changes of parameters. 



5.0.1 Estimating the Transition Probability for Resonant Mixing 

We turn to estimating the transition probability in the case of resonant mix- 
ing. We assume that the plasma frequency changes linearly along the path 
according to the relation Eq. [52j We again take the direction of the back- 
ground magnetic field fixed so that we consider the two state mixing problem. 
In this case we find, for example, 

d z b 2 = 6i- 29 * Bu e-^o tf(^-Ml) (68) 

° z2 1 2[{C-azf + {2g^Bujf} [ ™ } 

Now if we assume that the exponential is the most rapidly varying factor on 
the right hand side of the above equation then we can replace the coefficient 
by its average value along the path. We then find, 



x 



L 

-i r L 

e 



m;-wjj(0) ml-u 2 (L) 
tan tan 



„„2 _ , ,2fr\\ ™2 

b 2 (L) « 62(0) + i&i(0) T 



2g ( j > Buo 2g ( j ) Buj 



1 

[rnl-u>*(0))* + (2 g4) Bu) 



(69) 
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where we have used 
2g< t> Bu 



(m 



1 

La- 



tan 



2g <j> Bu 



— tan 



2g (j> Bu 



It is clear that the second term on the right hand side of Eq. [69] is small 
and hence the transition probability from one state to another is negligible 
as long as L is sufficiently large such that 



L » 
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1/4 - \r- 



(70) 



However this condition is not general and is applicable only in the limit of 
large uo. The basic problem is that close to the region where u p = m^, our 
original assumption that the coefficient of the exponential in the right hand 
side of Eq. 1681 varies slowly compared to the exponential term is not correct. 
We see this by comparing the magnitudes of the logarithmic derivatives of 
these two factors. We find 



D 1 = A log [{C 



az) 



-2a(C — az) 



d 2 = Aio ge -*/o'£(4-M a -) = jL 

dz 2u 
The assumption, \Di\ « \D 2 \, implies 



\2a(C — az) 



(C - az) 2 + (2 9( pBLuy 



(C - az) 2 + (2g^Bu) 



(C - az) 2 + (2g^Bu 
In regions where \C — az\ » 2g ( j ) Buj this condition demands 

lad << — (C - az) 2 . 
Since (C — az) is of order rr& this condition implies 



(71) 



\a\ « — mi . 

We next consider the limit where \C — az\ is small or comparable to 2g t f ) Buj. 
We test the condition, Eq. (TTJ at the point |C — az\ = f3(2g^Buj). In this 
case we find 

M « (/32+ ^ 1)3/2 (^) 2 ^ 
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This condition is most stringent for very small values of (3. However in 
this region, as we show below, the transition probability between instanta- 
neous eigenstates continues to be very small provided the condition Eq. [70] 
is obeyed. If we take (3 of order unity then this condition implies \a\ « 
(g^B) 2 uj, which can be expressed as, 

(72) 

which is rather stringent for small frequencies and magnetic fields. In re- 
gions where this condition is violated we can no longer ignore the transition 
between different instantaneous eigenstates. 

We can obtain an estimate of the transition probability in regions where 
Eq. [71] is violated. Assume that in these regions the exponential term in 
Eq. [68] varies much slowly compared to the coeffient. We then replace the 
exponential with its average over the path. We perform the integration over 
the region where C — az ranges from (3{2g^BuS) to —(3(2g^Buj). In this region 
we find 



Ab 2 = tan" 1 (3 &C4-m1) 



(73) 



The exponential term gives a contribution of order unity. It follows that the 
transition probability between different eigenstates is non- negligible, unless 
P«l. 



Summary The above results can summarized by stating that in the case 
of resonant mixing the adiabatic approximation is valid only if both the 
equations Eq. [70] and Eq. [72] are satisfied. If either of these equation is not 
respected then the transition probability between instantaneous eigenstates 
cannot be ignored. 



6 Applications 

The basic aim of the present paper is to make detailed predictions for polar- 
ization observables due to pseudoscalar photon mixing, which can be tested 
in future astrophysical and cosmological observations. Tests may involve 
CMBR, polarization observations from distant sources, or propagation of 
electromagnetic waves through regions of strong magnetic fields such as the 
pulsar magnetosphere or Active Galactic Nuclei. The resonance phenomenon 
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may play an important role in many of these situations since the waves prop- 
agate through regions with large variations in plasma density. A detailed 
study of all this phenomenon is beyond the scope of the present paper. Here 
we shall confine ourselves to simple estimates. 

We first determine the range of frequencies for which the adiabaticity 
condition is satisfied in the case of resonance in the supercluster magnetic 
fields. Using equation [721 with B = 0.1 fj,G [2D], n e = 10" 6 cm" 3 , L = 10 
Mpc and the current limit on the coupling g$ = 6 x 10 -11 GeV -1 , we find that 
adiabaticity is satisfied if uo » 5 x 10~ 3 eV. Hence we expect a significant 
effect at optical frequencies but negligible at radio. We should point out 
that if resonance occurs it gives a considerably enhancement of the mixing 
phenomenon in comparison to the case of uniform background. For the latter 
the mixing probability is limited by the oscillation length I in the medium and 
is proportional to {gBl) 2 [15]. For low frequencies this is very small and gives 
a significant effect only if u > 5 eV. On the other hand for a spatially varying 
medium if u p = somewhere along the path, the mixing probability is of 
order unity as long as the adiabaticity constraint Eq. [72] is satisfied. This 
constraint is satisfied at much smaller frequencies in comparison to what is 
required to obtain a significant effect if uo p ^ all along the path. 

For the galactic magnetic fields B pa 3 /iG and plasma density n e pa 
0.03 cm~ 3 and distance scale 50 Kpc, we find that adiabaticity condition is 
satisfied only for u » 50 eV. Hence here it is satisfied only for ultraviolet 
frequencies. Hence we can expect significant effects due to pseudoscalar- 
photon mixing at such high frequencies. 

The effects of resonance are most dramatic if (gBl) 2 << 1, where I is the 
oscillation length evaluated at some point far away from resonant region. In 
this case the resonant effects dominate. For supercluster magnetic fields this 
is found to be the case for a wide range of frequencies 5 x 10~ 3 < uj < 1 
eV (here the lower limit is determined by the adiabaticity condition). In 
this case the results Eq. [51] apply and the wave becomes almost entirely 
linearly polarized after crossing the resonant region, independent of its state 
of polarization at origin. If the direction of the magnetic field remains fixed 
during propagation then its state of polarization is well described by Eq. EU 
In general, however, the magnetic field will also twist along the path. In 
this case the wave will also acquire circular polarization during propagation. 
This phenomenon is discussed in a separate paper [29]. If the frequency is 
large, such that (gBl) 2 >> 1, then the effect of pseudoscalar-photon mixing 
is large irrespective of whether resonance occurs or not. 
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The pseudoscalar-photon mixing has been proposed [27] as a possible ex- 
planation of the observed large scale alignment of optical polarizations from 
distant quasars [30], [31] . The possibility of resonance in supercluster mag- 
netic fields further enhances the parameter space over which this explanation 
may be applicable. This explanation can be further tested by observing the 
spectral dependence of all the Stokes parameters. Pseudoscalar-photon mix- 
ing also predicts correlations between these parameters, which can be tested 
in future observations. 

The resonance phenomenon may also considerably enhance the produc- 
tion of pseudoscalars due to their mixing with photons as the electromagnetic 
wave propagates through magnetics fields in astrophysical objects such as the 
active galactic nuclei, pulsars, magnetars, galactic clusters etc. If we assume 
resonant production of pseudoscalars from active galactic nuclei at optical 
frequencies then these pseudoscalars can convert back into photons during 
propagation in the local supercluster. It is reasonable to assume that the 
direct photon flux, integrated over the entire source, would not be strongly 
polarized. Hence the observed polarization may be determined dominantly 
due to the conversion of pseudoscalars into photons in the local magnetic 
field. This may explain the observed alignment of optical polarizations from 
large redshift z > 1 quasars over the entire sky [31] as well as their correlation 
with the supercluster equatorial plane [30] . 



6.1 Cosmic Microwave Background Radiation 

We next consider the effect of supercluster magnetic field on the CMBR due 
to its mixing with pseudoscalars. Here even a small effect may be observable 
due to the high precision with which CMBR has already been measured. 

Intensity Effects: We first consider the intensity on the CMBR due to 
decay into pseudoscalars. Let P(9, 4>, E) be the probability of decay into 
pseudoscalars, where E is the energy of the photon and (9, 0) are the angu- 
lar coordinates. Let the original spectrum be denoted by f(E; T) and the 
distorted spectrum due to mixing with pseudoscalars by f'(E, T'). These are 
related by 

f(E;T') = f(E;T)(l-P(9,(j>,E)) (74) 
Assuming that P(9, <f), E) << 1 we find 

T'(9, 4>, E) « T(9, 4>) - ^'t^V - exp(-£/T )) (75) 
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where T is the mean CMBR temperature. As expected mixing with pseu- 
doscalars will produce both a frequency and angular dependence of the tem- 
perature. To be consistent with observations we expect P(6, <f), E) can have 
a maximum value of order 10 -5 . 

Angular Distributions: It is interesting to speculate that the mixing with 
pseudoscalars in the local supercluster might give rise to prefered orientation 
of the CMBR quadrupole and octupole [321 E3] due to a possible angular 
dependence of the mixing probability. This may also explain why higher 
multipoles are aligned with the CMBR dipole [31]. "Cosmic variance", in- 
voked to explain away magnitude puzzles, cannot credibly explain the mul- 
tiple coincidence of multipole directions. Hence an explanation in terms of 
local foreground effects is attractive. For consistency the mixing probability 
has to be of order 10~ 5 , which is roughly the strength of these multipole mo- 
ments. Current limits on the pseudoscalar photon coupling do not rule out 
this possibility. However if local effects are the source, why should they be of 
the same order, multipole by multipole, as cosmological CMBR fluctuations? 
Interaction with a coherent pseudoscalar field might be the explanation with 
the fewest number of arbitrary assumptions, although some fundamental re- 
vision of assumptions in cosmology might be needed to reconcile everything 
observed. Nevertheless it is interesting to investigate pseudoscalar-photon 
mixing in supercluster magnetic fields further since it may provide more 
stringent limits on the coupling g^. Pseudoscalar-photon mixing will also 
give rise to a spectral dependence to the CMBR temperature. Hence the ef- 
fect has to be sufficiently small so that it does not conflict with the observed 
agreement [35] of CMBR with black body radiation formula. 

Polarization: Pseudoscalar-photon mixing can also generate polarization 
of the CMBR radiation. If the mixing probability is of order 10~ 5 , which is 
required to explain the alignment of quadrupole and octupole, then it will 
generate polarization of the same order of magnitude. Hence if the polariza- 
tion is observed to be much smaller than this, it might be possible to impose 
further limits on the pseudoscalar photon coupling. The pseudoscalar-photon 
mixing can also change the CMBR polarization parameters without changing 
the overall degree of polarization. This effect is in general much larger. How- 
ever if the degree of polarization of CMBR is small, as expected, any change 
in the polarization will be observable by the current detectors only if it is of 
order unity. Hence here we require a relatively large contribution in order for 
new effects to be observable. The shift in polarization is of the order of the 
phase acquired by the wave after propagating a distance L [27]. For uniform 
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medium this is of order (g^B^lL, where I is the oscillation length. For the 
Virgo supercluster parameters and microwave frequencies of order 100 GHz 
we find that the phase is of order 0.1 for the coupling = 6 x 10~ n GeV -1 . 
This phase will generate all the Stokes parameters with relative strength of 
order 0.1, assuming the current limit on the coupling g^. Hence this contri- 
bution is small but it may be observable in future detectors. Similar results 
are found if we take into account the fluctuations in the plasma density [27] 
or propagation through intergalactic medium |17j . 

Condensate: A pseudoscalar condensate can also affect the CMB polar- 
ization [12]. In this case the rotation of polarization is equal to g<j,A(j), where 
A(f) is the total change of the pseudoscalar field along the trajectory of the 
electromagnetic wave. This effect, however, cannot generate circular polar- 
ization. This effect is independent of frequency and will also affect radio 
wave polarizations from distant sources [27J [31]. If this pseudoscalar field 
distribution is anisotropic, it will lead to an anisotropy in both the CMB and 
radio polarizations [36J. Future observations can provide stringent limits on 
this effect. 

7 Dark Energy 

Dark energy has come to denote a cosmic energy-momentum tensor of the 
vacuum. The source of dark energy is unknown. Unfortunately, the tradition- 
ally minimal option to employ one universal cosmological constant appears 
less and less credible. It appears more likely that dark energy is the grav- 
itational trace of a field, somehow chaperoning the evolution of Big Bang 
pressures, densities and phase transitions predicted by particle physics [37J. 
The need for a causal inflaton field is rather clearly manifested in the high 
uniformity of the cosmic microwave background (CMB) radiation. 

By now, cosmology as a whole perhaps cannot do without a dark energy- 
associated field. Yet of all interactions, exploring the Universe with gravity 
has a problem, in that gravity is the finest example where the coupling to 
fields is unknown! Unconventional as the remark may seem, in quantum 
theory there exists no way to find independently that part of the energy 
momentum tensor which couples to gravity. The rules of coupling are un- 
known, and if at all knowable, appear to hinge on ultra-high energy physics 
not experimentally testable. The method of using gravity as a probe, when 
it is not really known what gravity couples to, flies in the face of the long 
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successful tradition of using perturbatively stable interactions to study new- 
situations. It seems ironic that entire fields are built using gravity to explore 
the Universe's evolution when gravity is the least understood interaction. 

If we assume that dark energy is associated with a scalar field 0, such a 
field will also have a coupling to electromagnetism given in Eq. [TJ Hence 
this field might produce all the physical effects we discuss in this paper. 
Couplings of dark energy to ordinary matter fields as well as their limits 
are discussed by Carroll [19]. Yet the mass of such a conventional dark 
energy field is expected to be of order 10~ 33 eV, and then much smaller than 
the intergalactic plasma frequency. Assuming such values we do not expect 
the conditions for resonance to be applicable. One may, however, consider 
generalized models of dark energy. Indeed it turns out that by invoking a 
false vacuum it is possible to explain both dark matter and dark energy in 
terms of the invisible axion [381 EH] ■ in this case the mass of the background 
field can be much larger than currently assumed. Hence we cannot rule out 
the possibility that conditions for resonance may be applicable in some cases. 

8 Summary 

If one were to approach light and pseudoscalar field mixing "lightly," then 
there is a sequence of facile, dimensionally based arguments that could be 
used. First, since has dimensions of inverse mass, and electrodynamics 
is a theory with no scale, one might claim that all effects were relatively 
proportional to g^ou, and must vanish for g^uj « 1. This is false: the effects 
are cumulative, and observable for exceedingly small g^ « 10~ 12 GeV -1 . 
Next one might observe that with a background field B, the effects must be 
of relative size g^B/u, and vanish for large to. This is again false, as there 
are several other scales, sometimes u occurs in the numerator, and we find 
that for realistic parameters there are observable effects persisting all the 
way from radio to optical frequencies. 

Then we return to our original goal of probing pseudoscalar field with 
light. It is certainly very interesting that light can be dimmed by interac- 
tions with pseudoscalars. Yet in comparison the variety and variability of 
polarization-based observables seems almost unlimited. We believe that as 
polarization observations accumulate, more and more anomalies will appear. 
It would be gratifying to have data so that a new and systematic study of 
pseudoscalars via a coupling that is stable under perturbation theory can 
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commence. 
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